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ABSTRACT 


The major environmental Ultra Low Frequency (ULF), magnetic 
noise sources observedon the airborne AN/ASQ-81(V) optically 
pumped magnetometer are electromagnetic micrepulsatrans, geo- 
logical structures and ocean wave noise. 

Experimental determinations of the amount of magnetic noise 
produced by ocean surface gravity waves were made. The author's 
current research designed to remove the wave noise and increase 


the probability of Magnetic Anomaly Detection (MAD) is summarized. 
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I. INTRODUCTION 


Magnetometry has been used as a classical geophysical pro- 
specting tool for years (Dobrin, 1952). The principles of 
magnetic anamoly analysis has also been applied to Anti- 
Submarine Warfare (ASW) with a limited amount of success for 
submarine localization. 

The sensitivity of magnetometers has been continually 
improving so that it is now possible ta observe, from an air- 


craft, fluctuatious of 107? 


nano Tesla/hertz (nT/Hz), imposed 
on the earth's constant field of approximately 50,000 nT. 
Land stationed superconducting magnetic snesors have demon- 


“6 it/Hz. The advent of airborne 


strated a sensitivity of 10 
superconducting magnetometers would greatly increase detection 
ranges, however, the extent of improvement depends in detail 
on the amplitude and frequency of the environmental noise which 
is above the noise floor of the instrument. 

Airborne Magnetic Anomaly Detection (MAD) for ASW is done 
in the .04 to 2 Hz frequency band. Ocean surface gravity waves, 
in the .03 to 10 Hz frequency band (Silvester, 1974) are con- 
tributors to the magnetic noise spectrum. The other two major 
environmental noise sources are electromagnetic pulsations and 
geological noise. 

It is generally assumed that noise due to ocean gravity 
waves is not "loud" enough to be of any consequence. However, 


experimentally it was observed that noise levels induced by 
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typical sea condition in the mid latitude oceans during the 
late fall to early spring are at least one order magnitude 

greater than the noise floor of airborne magnetometers. The 
noise greatly reduces the effectiveness of Magnetic Anomaly 


recognition. 
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II. THEORY 


The problem of determining the electric current induced 
by deep water gravity waves at the ocean surface and the resul- 
tant periodic fluctuations in the magnetic field has been de- 
scribed by Weaver [1965]. The magnetic field fluctuations 
above and below the ocean were obtained by solving Maxwell's 
Equations for the magnetic flux density vector in air and water, 
and applying electromagnetic boundary conditions at the sea 
surface. 

The movement of the sea surface through the earth's mag- 
netic field induces electric eddy currents in the sea water. 
In deep water, assuming linear ocean wave theory, the water 
particle motion is periodic and circular, as is the induced 
Magnetic field. This process becomes more complicated as one 
decomposes the sea surface spectrum into waves of various 
frequencies which may be propagating in varicus directions 
with respect to the magnetic field vector. Further comrlica- 
cions arises as nigher order wave theories are applied to de- 
scribe the water particle motion and the waves propagate into 


shallow water. 


A. MATHEMATICAL MODEL 
The mathematical model of the induced field was developed 
by Weaver [1965], assuming linear wave thecry to describe the 


water particle motion. 


It was derived for a monochromatic ocean surface gravity 


wave, oropagating in deep water, in the X direction. The 
development is summarized below. 
The coordinate system formed a right hand triad as follows: 
X-axis, was taken in the direction of wave propagation. 
Y-axis, was taken parallel to the wave crest. 
Z-axis, was taken. downward into the ae Z=0 at the 
still water level (SWL). | 
It was assumed that the earth's magnetic field was constant 
in the region, and the wave crest was infinitely long. The 
model tecame two dimensional as all magnetic gradients in Y% 
were negligible. 
Z The earth's vector magnetic bield, Bi, was expressed as 


mn 


| eee, rN 
B, = |B. [feos cost - cesT smi Oj + sin k) 


aA a 


r 


where i, 47, k were unit vectors alcng the three coordinate 
axis x, y, Zz. I, was the magnetic dip angle in the vertical 
dlane, relative to the horizontal. @, was the inclination 
in the horizontal plane of the direction of wave prcepagation, 
clockwise from the magnetic meridian (see Figure 1). 
The field vectors satisfied Maxwell's equation in the form 
= R 
Wx = - 2B ) (1) 
“4 
— a : xs 
Wx8 = cEt+ 6(Vx(B, +B) + es. 
Me dc 


(2) 


XA. 
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= . y -_ 
where © was the electric field vector ( /Meter), B was the 


induced magnetic flux density (Tesla), H, was the permeability 


of free space (unx lon’ Henry Meter), Oo was the electrical con- 


(Siemen 


a 
ductivity /Meter), V was the water particle velocity 


vector (Meter 


/sec), € was the permittivity (Farads peter), and 
t was time (seconds). 

The deep water particle velocity, described by linear wave 
theory, was of the form 
ie va. He 6 (We -mx) “my 
VeAu(ic+R)2 


where a was the wave amplitude (Meters), w=2n/f was the angular 


Radians 


frequency of the wave ( /sec), mz2n/L was the wave number 


(Radians py eer) , and t was time (seconds). 
—_ 
Since V was harmonic in x and t, the magnetic and electrical 


=> 
field resulting from V has a similar structure, thus 


~ = L (we =m) 
Be bye ; (4) 
=- = L(wt -mx) 
E=*Emz : (5) 


The 8 field and the F field are nonedivergent which implied 


“ObG) 
= ae Dy (4) . (6) 


g 


Spit = Lm Eng . (7) | | 


The differential equation in water satisfied by B was 
obtained by assuming B<<B_ and we<<o. Then, taking the curl 
of equation (2) and using the non-divergent portion of the 
vector identity, Ux(0xk) = -7°A + 9(0:A), the following was 


obtained: 


~V7B eOUXE+ CHVX(VXB) . (8) 


Further simplification of equation (8) was obtained by the 
elimination of VxE by equation (1). From equation (3), V was 
eliminated. Then if y? = iu wo, the following equation results, 

*h 2 

aby 2 2 2 . . -™, 

2 = (g'rm?) ba -Jiam (Buon cBuyNit+R)e a (9) 

In region Z>0, equation (9) was solved for the Z component 
by assuming 


dato) =Qe (rent + am(B«)- iBuwle "7 ; 


where Q was an undetermined coefficient. 
In air €>>0, thus equation (2) became, 
B. ook 
o 
The curl of euation (10), when FE was eliminated by equation 


(1), resulted in 


pis) a( m= 4&0) bey ; 


D3" (11) 


Confining the sea surface to scales of motion described by 


gravity waves implied that m>>y ew’, thus 


13 


oP SLT SRST Oe ee BET a LO ee a, Be ee ee a a eT ee oe a a ee ee ee a SE Dy ee SEE OE a ee ERT Le a eT ee 


~—> 


a —> 
2 bu» 2 mtb, ; (12) 
ox 


In the region Z<0 the solution for the Z component of 


equation (12) was of the form 
fa (13) 


where R was an undetermined coefficient. 
The magnetic field vector was contint.ous across the air-sea 
interface, thus, the following boundary condition at Z=0, existed: 


Dy aa (20) = b, waren, (2.20) 


Ob, ae (2 20) - ‘Ob; wan (220) 
oy 4 (14.d) 


(14.a) 


Applying the above boundary values to equations (9) and (13) 


yielded R= q rf am (Bun = { Pet>) ) 


(15) 


and 
mR = -G(se nt) = am? (Ban -i Bip). (1s) 
Equations (15) and (16) lead to 
Qe rZmta( Pein - (Beg) 
ym +(e 


(17) 
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and 


R= % (Bem) =e Bary yemt ym) 
es 
x 


Thus the Z component of the induced magnetic field vector 


(18) 


in air was Zz 
Dof3) > Qm ( Boo +e Boy \Cus Pers) 2 a 
' en nanan na : 
os 
f2 


In order to further simplify equation (19), << was true 
m 


(19) 


for gravity waves resulting in 


zu 2 
_ af : (20) 
by ” Tn (Pron -(Bege3 ; 
From equation (6) and (20) it was seen that 
COxcy) = bey) (21) 
Tt was convenient to measure the earth's field with a 


total field magnetometer. Thus, 
Besensveny 7 Bet 6,69) cosleesErbyey sw, (22) 


when equations (20) and (21) were used to simplify equation 
(22), and only the perturbations were considered, it can be 
shown that the magnitude of the observed magnetic field at a 


specific frequency would be, 


b(F) = AGH | (cosz ese) +(su ry) e"9 4 (23) 
ZL0, 
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The deep water wave length (L) of a surface wave is related 


to the wave period (T) by, 
z 
{= qt (24) 
ya 
where g is the acceleration due to gravity. 


Eliminating m, from equation (23) by the use of equation (24) 


leads to, 


tf) = seth [Rl(osz cos oy +(sin Dp 37 g 8) 


Equation (25) was used to model the magnetic field spectrum 
when the directional ocean wave spectrum was known. 
Using an airborne magnetometer to measure the magnetic 


Spectra, there was a doppler type frequency shift, described as, 


f’ = #(i- (2m fu/g) os ¥) ae 


(see Appendix C). Where f' was the wave frequency observed on 


the aircraft (sec}), £ was the actual wave frequency, u was 


(Meter 


the velocity of the aircraft /sec), and w was the angle 


between aircraft motion and direction of wave propagation (see 


Figure 1). 


i 
‘ 


III. THE EXPERIMENT 


Magnetic data was collected by an airborne magnetometer 
flown over the sea surface in order to verify the theory pre- 
viously presented and to further understand the observed environ- 
ment in the 0 to 2Hz band. The AN/ASQ-81 magnetometer, aboard 
the P=3C Orion aircraft, is capable of measuring the signal 
produced by ocean gravity waves which are frequency shifted 
due to relative motion between the wave and the aircraft flying 
at 100 meters altitude. The analog voltage time series output 
of the magnetometer was recorded on an HP3960, FM tape recorder 
for post flight analysis. 

The fluctuation of the sea surface below the aircraft was 
determined by two techniques. The analog voltage time series 
output of the APN-194 radar altimeter, also aboard the aircraft, 
was recorded. 

Independent of the aircraft data the ocean surface wave 
climate was determined from analysis of the Fleet Numerical 
Weather Central (FNWC), spectral Ocean Wave Model (SOWM) and 
Measured wave spectrum from National Oceanic and Atmospheric 
Administration (NOAA), Environmental Buoy, NR-6003 (EB-17), 
located at 52.2N, 155.8W from the airborne data. { 

The airborne data was collected in vicinity of 52°N, 156°", 


around the buoys' reporting period, 09Z, 10 Feb 78. 


A. AIRBORNE DATA COLLECTION 
The AN/ASQ~81 employs the atomic properties of optically 


pumped, metastable helium atoms to detect variations in total 


17 } 
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Figure 2. Region of data collection. (xX) FNWC Grid Points, 
(&4)NOAA Buoy, (---) First Leg of Aircraft Data, 
(+--+) Second Leg of Aircraft Data. 


18 


jira hd Bre os 


field magnetic amplitudes. The output of the sensing element 
is converted to an analog voltage which is bandpassed and scaled 


to display field variations (Order 10972 


T) as a function of time. 
Instrument specifications are described in Appendix A. 


Perpendicular legs, 270° True and 000° True, were flown at 


"Gach ¢ 


approximately 60 meters altitude and 180 knots (91.44% M/s) for 


(i 


30 minutes each. (See Figure 2) The magnetic data was collected 
in the .04 Hz to 2 Hz band, thus, eliminating the earth's con- 
stant field and any micropulsations above 2 Hz. Aircraft motions 
were compensated for by the magnetometer using baroaltimetry, 
selective filtering and accelerometers. No variations in the 
magnetic field due to the aircraft were expected due to the 

above real time processing and the straight and level flight 
paths flown. 

The instantaneous aircraft altitudes were recorded simul- 
taneously with the magnetic data from the APN-134 radar altimeter, 
which responded linearly to apparent altitude changes due to 
long period surface gravity waves. The wave frequency was 
shifted due to aircraft motion (see equation (26)). Beam 
characteristics, 60° ecne, were such that at 60 meters the 
diameter of the radar spot on the surface was approximately 70 
meters. In order to resolve wave amplitudes from 60 meter 
altitude, the wave period must be greater than 10 seconds. 

The frequency and amplitude errors attributable to the size 
of the radar spot were considered in the analysis. 

The AN/ASQ=81 and APN-194 signal pick-off points are 


described in Appendix A. 


3 
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B. SPECTRAL WAVE DATA COLLECTION 

The experiment was conducted around synoptic reporting 
period 00Z, 10 Feb 1978 in vicinity of NOAA Buoy, E2-17, so 
the spectral information computed by the Buoy for the 002 time 
period was useful in interpreting the aircraft data. 

In addition to Buoy data the SOWM output was used in order 
to infer directionality and improve the spectral wave estimate 
in the region of the flight. 

The Buoy is designed to report meteroslogical and oceano- 
graphic conditions at the sea surface every 3 hours. The wave 
conditions are measured for 20 minutes by accelerometers and 
reported for the Buoy location as frequency and amplitude 
(Steele and Johnson, 1977). The Buoy data of 902, 10 Feb 1978 


is reported in Table I. 


Frequency (Hz) Period (sec) Amplitude (meters) 

»06 16.7 8322 

G7 14.3 27574 

08 12.5 1.6064 

09 11.1 1.9604 

pat 10.0 1.9312 

«ll 9.1 1.5419 

+12 8.3 1.0116 

-13 Ted -9083 

-14 7.1 .7917 

15 6.7 6543 | 
.16 6.3 4793 . 
+17 5.9 -4$11 

.18 5.6 4212 : 

Table I 


NOAA Data Buoy 002, 10 Feb 1978 
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The FNWC, SOWM uses the analyzed surface pressure field 
to determine wind velocities and then builds, radiates and dis- 
sipates waves as prescribed by wave theory (FNWC Tech Note No. 
75-3). Model output is energy as a function of frequency 
and direction at each grid point. The wave energy (E) at a 
frequency can be converted to wave amplitude (a) by a= E-df. 
SOWM frid coverage of the data area is shown in Figure 2. 


The model output for the closest grid point is summarized in 


Table II. 
DIRECTION FROM 283° 313° 343° 013° 
(OT) 
PERIOD (SEC) AMPLITUDES (METERS) 
16.4 .347 ) .598 £987 
14.93 .578 .408 1.902 £578 
13.89 .708 1.293 1.082 £347 
1233 1.636 0 1.686 1.686 
12.04 1.027 1.686 sie .959 
10.99 Leu? 1.584 314 74192 
teal 1.993 1.493 944 LLu36 
8.5 £944 1.356 1.227 ee we 
7:5 6851 £973 987 .914 
6.5 2409 2457 6354 167 
Table II 


FNWC, SOWM Data, 032 10Feb78. Grid point located at 
$2.357N, 154.099W. Amplitudes (meters) 
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IV. DATA ANALYSIS 


A. TIME SERIES ANALYSIS 

The general form of the magnetic and altitude time series 
is shown in Figure 3, which is a sample of each time series. 

The time series were digitized at .1 second intervals, 
resulting in a Nyquist frequency of S Hz. Any spectral power 
above 5 Hz would be folded back into the spectral estimate 
(alliasing). The magnetic data, due to the low pass, 2 Hz, 
characteristics of the instrument, alliased only instrument 
noise which was at a low level. It was assumed that no sig- 
nificant alliasing of the altitude data occurred due to the 
altimeters 0.1 second RC constant which results in a 1.59 Hz 
cutoff frequency Cw, = 1iRC). 

Prior to computing power spectral estimates, the data was 
calibrated to insure constancy of unit magnitudes. The magnetic 
data voltages were linearly calibrated to nT and the altimetry 
voltages were linearly calibrated to meters after the flight. 

The data was prepared for the Fast Fourier Transform (FFT), 
(Bendat and Piersol, 1971) by forming 16 ensembles, 192.4 seconds 
each. Each ensemble was then linearly detrended. The detrending 
removed significant variance from the frequencies below .001 Hz. 
The cosine function was used to taper 10% of the time series 
by multiplying the first and last 51 points of the time series 
by the cosine. Points 1 and 1024 become zero and points 52 
thru 973 are as digitized. The tapering procedure is done in 


order to prevent non-real low frequency power leakage, due to 
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Horizontal 


Vertical scale is lmm/volt. 


Sample of magnetic time series (top) and 


altitude time series (bottom). 


scale is lmm/sec. 


Figure 3. 


ensemble averaging, into the spectrum. The variance removed 
in the time series taper is restored to the power spectra after 
the FFT (Bendat and Piersol, 1971). 

The FFT computes the real and imaginary components of the 
Fourier Transform. It was then possible to compute the sample 
power spectra of each signal by ensemble averaging. The normal- 
ized cross power spectra, defined as coherence and the phase 
of the cross power spectrum between signals were also computed. 

Each of the two flight courses were treated as separate 
data sets. Note in Figure 4 the power spectra of the altimeter 
(4? siz) and magnetometer (nt? Hz) and the coherence and phase 
as a function of frequency. The numerical results are summarized 


in Appendix D. 


-B. SPECTRAL WAVE ANALYSIS 

The analysis of the rNWC, SOWM directional wave spectra 
and EB-l7 ocean wave amplitude spectra in order to describe 
the magnetic field variations observed aboard the aircraft, 
required fitting the large scala information to the microscale. 
This change of scale was feasible due to spatial uniformity 
of the ocean wave spectra over the area and frequency range 
investigated. 

The ocean wave climate in an area can be described by that 
at a point over a period of time. £B-17 reported the waves 
during a 20 minute period to have the spectral shape listed in 
Table I. The FNWC ccean wave amplitude and directional spectra 


had a Similar shape listed in Table II. 
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The directionality of wave components at E£B-17 was deter- 
mined by group velocity analysis (Silvester, 1974) which pro- 
pagated the FNWC, SCWM data at nearby grid points to the buoy 
position. The greatest weight was On the grid point closest 
to the position of the buoy. The grid point was located at 
52.357N, 154.099W, 60 miles from the buoy. The analyzed buoy 


directional amplitude spectra is reported in Table III. 


PERIOD (SEC) AMPLITUDE (METERS) DIRECTION FROM (°T) 
1657 .8322 340° 
14.3 .7574 340° 
2.05 1.6964 315° 
S140 1.9604 315° 
10.0 1.9312 320° 
9.2 1.5419 420° 
8.3 1.0116 320° . 

ee 9083 320° 
74 £7917 330° 
ee 6543 330° 
6.3 £4793 320° 
5.3 24511 320° 
5.6 212 320° 

& 

Table iI! 


Analyzed directional wave data at NOAA buoy E8-17, 00Z 10 Feb 73. 


It is important to note that the NOAA data and the FNWC 
data have similar spectral shapes although amplitudes are not 
equal, verifying scale analysis assumptions. The buoy reports | 
a spectral peak at approximately 11 seconds. FNWC predicts a 


Spectral peak at approximately 12 seconds. 
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The observed ocean wave climate on the aircraft, although 
similar, was not the climate observed by EB-17 or predicted 
by SOWM. Each of the ocean wave data sources had errors asso- 
ciated with them, which would contribute to the differences 
in the data. Scale considerations, length of airborne data 
record, modeled grid point locations, buoy location and flight 
path, further contribute to the uniqueness of the separate 
data sources. 

Due to the multidirectional nature of the SOWM data and 
its proximity to the flight path, it was assumed that the 
closest FNWC data point (Table II) most closely resembled the 
ocean wave climate observed on the aircraft. The SOWM wave 
information was applied to equation (25) resulting in an ex- 
pected magnetic amplitude spectrum. The eypected frequency 
shift due to the aircrafts motion relative to the waves was 
computed by applying equation (26). To be consistent with the 
power spectral density computed from the aircraft data the 
amplitude spectrum described above was squared and normalized 
by multiplying by the frequency interval of each frequency 
estimate. The results are summarized in Figure 5 and numeri- 


cally in Appendix D. 
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V. DISCUSSION 


Baker and Graefe (1968) showed experimentally that the 
theory reviewed in Section II describes the magnetic field 
induced by ocean surface gravity waves. Comparison of the 
aircraft magnetic power spectra, Figure 4, and the expected 
magnetic power spectra computed from SOWM data, Figure 5, 
further verifies the electromagnetic theory described above. 

Figure 6, which is the superpositioning of Figure 4a on 
Figure Sa and Figure 4b on- Figure Sb, illustrates that although 
the aircraft and wave data magnetic spectra are not the same 
they are similar in amplitude and frequency. The differences 
observed in the low frequency end can, in part, be attributed 
to the assumptions nedessary to the formulation of each magnetic 
power spectrum and the assumption that each data source is 
stationary in time and space. 

The FNWC data, although close to the flight path, was not 
coincident with the aircraft track (see Figure 2). Furthermore, 
determining the barometeric surface pressure field over data 
Sparse areas, such as the Gulf of Alaska caused inaccuracies 
in the SOWM wave height, frequency and directional information, 
which could have contributed to some of the differences observed. 

When studying the magnetic spectrum computed from the air- 
craft data it cannot be assumed that the only source of noise 
was the ocean gravity waves. Although no significant electro- 
magnetic micropulsations in the 9 to 2 Hz frequency band, were 


expected, they cannot be ignored. 
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Production of magnetic noise by the aircraft's motion and 
onboard equipage was considered. Ideally, most aircraft pro- 
duced noises were compensated for in the magnetometers real 
time processing scheme and the flight was straight and level 
on each leg which would further reduce noise. However, expe- 
rience shows that the noise sources described above are not 
completely remcved from the spectral analysis when long data 
records are required. 

Note in Figure 6a that both the aircraft and the analyzed 
SOWM magnetic power spectra indicate a broad coincident frequency 
band. In Figure 6b it is seen that a change in the flight path 
causes the band width of coincidence to compress, shift to 
higher frequency, and two distinet maxima are observed in the 
aircraft spectrum. One of which is nearly coincident with the 
SOWM magnetic spectrum. This indicates that a significant mag- 
Netic noise source is external «5 the aircraft and ccean surface 
wave dependent. 

Most significant, is the infcrmation implied from the analysis 
of Figure 4 and Figure 6. It is seen (Figure 4) that the co- 
herence is high and the phase relationship is good between the 
magnetometer and the altimeter in a small frequency band (.2 to 
-7 Hz). The trend observed in the phase relationship over that 
band can be attributed to signal delays imposed by the electrical 
characteristics of each instrument. 

This highly correlated frequency band falls in the same 


band in which the ocean gravity wave noise is expected. The 
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altimeter specifications indicated that it was accurate, + l 

foot and has a linear response to about 1.5 Hz. Thus it can 

be concluded that the altimeter senses the sea surface flucuations. 
Note the large spectral peak below .l1 hertz on the alti- 

meter power spectra of Figure 4. The peak is approximately 

one order of magnitude greater in the data of the first flight 

leg (Figure 4a) when compared to the seccnd (Figure 4b). This 

peak was due to an intermittent square wave caused by malfunc- 

tion of the altimeter. The square wave does not significantly 


degrade the analysis presented above. 


VI. APPLICATION 


It was shown in section V that the altimeter can be used 
to describe the instantaneous fluctuations of the sea surface 
below the aircraft which are directly related to a portion of 
the magnetic noise. Thus, the high correlation between the 
altimeter and the magretometer indicates that it is possible 
to remove a portion of the magnetic noise due to ocean surface 
gravity waves. 

Current Naval ASW tactics require that the magnetic field 
information be analyzed in the real time domain. Thus, noise 
cancelling must be done real time. 

Spectral analysis or convolution filtering is not practical 
due to the signal characteristics of the target anomaly and 
the long length of records required to compute the FFT which 
describes the low frequency signals. 

The magnetic anomaly is a single transient signal imposed 
on the total magnetic field thus making spectral analysis 
impossible. Furthermore, it is possible that the aircraft 
may change direction and flight attitude before a sufficiently 
long record is collected from the altimeter and the magneto- 
meter to determine the coefficients of the Fourier Transform 
for convolution filtering. 

In order to keep noise removal in the real time domain the 
principle of Adaptive Noise Cancelling (Widrow, et.al., 1975) 


was applied to the 10 Feb. data. 


The adaptive noise cancelling was done by a microprocessor 
developed at Naval Oceans System Center (NOSC) and already 
used in frequency analysis at frequencies higher than those 
analyzed for ultra low frequency magnetic information (0-2 Hz). 
The altimeter (reference signal) contained noise that was corre- 
lated with noise in the magnetometer (primary signal). Each 
signal was digitized and the algorithm used computed the solu- 
tion which removed the correlated noise from the primary by 
adaptively subtracting the reference signal to converge on 
minimum power in the primary signal. 

Preliminary results indicate that this method mey be useful 
for real time noise cancelling in airborne magnetic data (see 
Figure 7 and Figure 8) while the magnetic anomalies not due 
to wave motion snould be relatively unaffected. 

Note in Figure 7 the relative signal levels before and 
after noise cancelling. The signal appears to be cut in gain 
approximately four to one after noise cancelling. However, 
spectral analysis (Figure 8) indicates that the only power 
reduction was in the frequency band that has been indentified 


as the ocean surface gravity wave noise frequency band. 
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Figure 8. Aircraft magnetic spectrum (A) compared to "noise 
cancelled" magnetic spectrum (3). When altimeter 
was used as noise reference. Zach spectrum has 
16 degrees of freedom. Note the band limited 
noise reduction indicated by the cross hatched 
area. 
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VIT. SUMMARY 


It has been Shown that ocean surface gravity waves con- 
tribute significantly to the environmental background noise 
observed in airborne magnetometry. The noise is the result 
of the magnetic field associated with electrical eddy currents 
induced by ocean wave propagation through the earth's magnetic 
field. 

The APN-194 altimeter aboard the P-3C ORION aircraft is 
capable of describing the wave climate below the aircraft sim- 
ultaneously with total field magnetic measurements observed 
by the AN/ASQ-81(V) magnetometer, also aboard the P-3C. Thus, 
the altimeter can be used to infer the wave noise measured by 
the magnetometer. 

Analysis demonstrated that it is possible to remove a sig- 
nificant portion of the magnetic noise due to ocean gravity 
waves by using the altimeter as a noise reference and sub- 
tracting that processed reference in real time from the mag- 
netometer time series output. 

Further research with real submarine targets is underway 
in order that the ASW community will benefit from applications 
of this work. Development of a greatly improved magretometry 


system is possible in the near future. 
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APPENDIX A 
MAGNETOMETER AND ALTIMETER SYSTEM CHARACTERISTICS 


I. MAGNETOMETER SYSTEM PARAMETERS - 
AN/ASW-81(V) MAGNETOMETER GENERAL CHARACTERISTICS 


Performance-Magnetometer Group 


Range 23,000 to 75,000 gammas 
Noise level 0.016 gamma (peak-to-peak) 
Figure of merit (FOM) 0.24 gamma 


MAD bandpass (selectable) 


High Pass Corner frequencies 9.04, 0.96, 0.08, 
and G.1 Hz 
Rejection slope 35 ¢B/octave 


Low Pass Corner frequencies 6.2, 0.4, 0.6, 
and 2.0 Hz 
Rejection slope 24 dB/octave 


Gain selection 0.1, 0.2, 0.4, 1, 2, 4, 10, 20, 
(full scale) and 40 gammas 
Altitude compensation Barometric type 


Reference altitude, self-zeroing 
Linear range, + 200 feet 


Environmental Limits 


tion MIL-=E=5400 within the limits described below. 
Temperature 


| 
This equipment meets the requirements of military specifica- 
! (1) Equipment Operating: -S4°C to 71°C except the 
Magnetic Detector DT=-323/ASQ-81(V) or Towed Body 
TB-623/ASQ-81(V) which have +55°C maximum 
temperature. 
(2) Equipment Nonoperating: -62°C to +85°C except the 
Magnetic Detector DT-323/A5Q-81¢V) or Towed Body 
TB-623/ASQ-81(V) which have +71°C maximum 
temperature. 


Altitude Sea level to 30,200 feet 


Humidity 0 to 100 percent 


Shack Up to 18 g with duration of 11 
+1 milliseconds 


eyrve IV (10 g, § to 500 H2) except 
Detecting Set Control C-6983/ASQ-81(V) 
or CH9086/ASQ-81(V) which are curve 
III (5 g, 5 to $00 Hz). 


Vibration 


Input Power Requirements 
Magnetometer Group 
115/200 volts ac 3-phase power 
Voltage limits 108 to 118 volts 
Frequency limits 389 to 420 Hz 


Power 150 VA maximum 


Panel Lamp Power 


Voltage limit 0 to 28.5 volts 
Frequency de or ac 
Power 10 watts maximum 


Mangetometer Pick off Points 


Signal pick off soints in the P-3C were Rack C-+3, T8-433, 


pins C-14 and C-15. 


Pa 


a eae 


II. ALTIMETER SYSTEM PARAMETERS - 
AN/APN-194 RADAR ALTIMETER 


Manufactured by Honeywell, Inc., Government and Ared Products 
Division 


Specifications 
#1 ft. accuracy to 40,000 ft 
60° beam width to 3d4B point 
-l sec RC constant 
25-35 nano second pulse 
0-25 nsec gate 
Altimeter Pick off Points 
Signal pick off in the P-3C were FS 1005, TB-416, pins 


A-5 and A-6. 
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APPENDIX B 
SPECTRAL WAVE CLIMATE 
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APPENDIX C 
DERIVATION OF SURFACE WAVE FREQUENCY SHIFT DUE 
TO MOTION OF AIRCRAFT RELATIVE TO THE WAVE 


The wave equation due to motion of the aircraft is, 


iwtemx~* iw “tem 
erutomx > ei x 


Where w is the angular frequency of the wave, mis the 
wave number, t is the time and x is the distance from a re- 
frence point in the direction of wave propagation. The primed 


terms are the values due to motion, 
x* = xt(Cucosy)t 


Where u is the aircraft velocity andyis the angle clock- 


wise from the horizontal wave motion vector to the horizontal 


t 


aircraft motion vector. 


Then, 


a iwt-m(x+Cucosy) t) é giv “i-mx 


Mell Sass 


which leads to 


2it-mul(cosyp)t 


then solving for wo” leads to 
w * _ = Wemu cosy 
when w = 2m re m = an: and L = al oa: 
equation 25 results 


£* a £.[1-< amfu/ cosy] (26) 
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; 
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APPENDIX D 
SPECTRAL VALUES OF AIRCRAFT DATA AND MODELED NOISE 


I. MAGNETIC CONSTANTS AT S2N, 156N FOR FEBRUARY 1978 


Earth's total magnetic field, B, = $0,807 nT 
° 


Dip angle, I = 20 


¥ ned Bites ete, atone 9.2 
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THE MODELED MAGNETIC NOISE PRODUCED FROM 
THE PREDICTED SOWM SPECTRA DATA LEG 1 


Table Del.a 


The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircraft Course 
270-. 


FREQUENCY (Hz) AMPLITUDE CT siz) 
5042 0002 
£03 0002 
£12 0024 
14 .0037 
SiG 0023 
.18 .99uL 
.20 0014 
2 .0001 
123 00064 
26 0015 
27 0028 
£29 .0102 
Pace 0039 
.36 .0042 
a7 .0076 ' 
142 0128 
4S 0048 : 
47 0054 i 
i 56 .9052 : 
7 0016 : 
i £91 .0018 
Ta. 00002 


FREQUENCY (H2) 


Biel 
13 
14 
.16 
.18 
21 
£25 
£27 
£28 
mee 
.32 
36 
37 
Ay: 
43 


1.13 


2 
AMPLITUDE ("7 sz) 


.00087 
.0012 
.0012 
.0037 
.0915 
-0014 
.0013 
.0057 
.0088 
00004 
.0087 
.0069 
.0063 
.60339 
.0096 
.0034 
0012 
0022 
~0002 
.00057 
.00919 
-00018 
.009013 


Table Del.b The Modeled Magnetic Noise Produced from 
the Predicted SOWM Data, Aircarft Course 
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Figure D-1 Altitude Spectrum (CS 000°). 


°3,00 


Le as‘t o8°6 et"o ore oe Ge°e- ae 


zH/ ALISNAG IWALIAdS 901 
z 


i 4 
J 
& 
se 
& 
¢ 
N 
x= 
~ 
GN 
& 
: : 
. ng 
fe 
ray 
oe 
Fx 2 
a . 
‘ z 
g i 
5 
m3 8 
A. < 
i) ‘ ms 
ro) : 
2 ? 
a 
8 i 
e t 
‘ ? 
: 
8 4 
Af j 
7 2.20 =1.60 ~0.90 30.20 0.50 F 


LOG FREQUENCY (Hz? 


Figure D-2 Magnetic Spectrum (CS 000°). 
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Figure De3 Altitude Spectrum. 
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Figure D-4+ Magnetic Spectrum. 
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